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Abstract 

A global-in-time existence theorem for classical solutions of the Vlasow-Darwin 
system is given under the assumption of smallness of the initial data. Furthermore it 
is shown that in case of spherical symmetry the system degenerates to the relativistic 
Vlasov-Poisson system. 

1 Introduction 

Kinetic models play an increasingly important role in todays plasma physics. On the one 
hand much effort is used to deepen our analytical understanding of some problems where 
no other description seems to be adequate. On the other hand progress has also been 
achieved especially with numerical simulations (see e.g. 

In the kinetic picture, the particle distribution of a one-species plasma is described by 
a time dependent density function f(t,x,p) on phase space. If collisions of the particles are 
neglected and a relativistic model is used, then / is subject to the transport equation 

d t f + v(p)V x f + K(t,x)X7 p f = (1.1) 

with force term K = E + v x B. Here E and B denote the electric and the magnetic field 
respectively and the relativistic velocity is given by 

v(p)= j=L= p. (1.2) 

Note that all physical constants such as the speed of light or the rest mass of the particles 
have been set equal to unity. 
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Eq. (jl.ip is usually called the Vlasov equation. Expressions for the charge and current 
densities p and j in terms of the phase space density / are given by 



To obtain a self consistent closed system one has to take into account how the ensemble 
modeled by the density / creates the fields E and B. Usually this is done with the full 
system of Maxwell's equations, but numerical difficulties of simulations of that system have 
stimulated a search for alternatives (compare [3])- The present paper deals with what is 
known as the the Darwin approximation. Here the electric field is split into a transverse 
and a longitudinal component as follows: 



In the evolution part of the Maxwell equations the transverse part of the electric field is 
neglected, resulting in 



The system consisting of the equations (11. ip - (11, 6h is called the Vlasov- Darwin system. 
The main feature of this system is that the field equations are elliptic which in particular 
facilitates a numerical treatment since a time integration step, which is needed to solve the 
Maxwell system can be avoided here (|llj). The justification of the model seems possible 
in case the particle velocities are not too fast when compared to the speed of light. 

Up to now there are only few mathematical results known for this system. In 2003 
Benachour et al. [2] proved an existence theorem for small initial data: this assumption 
implies global-in-time existence of weak solutions of the Cauchy problem. Later Pallard 
[8] removed the smallness assumption and added a result about solvability of the Cauchy 
problem in a classical sense: To a given initial datum /o £ C|(M 6 ) there exists T > and a 
classical solution /: [0,T[xR 6 -^R of the Vlasov-Darwin systems satisfying f(0) = Jo- 
in the main part of the paper we present a result which is well known for the Vlasov- 
Poisson system (VP), the Relativistic Vlasov-Maxwell system (RVM), and other related 
systems such as the spherically symmetric Vlasov-Einstein system (cf. [HE! [10]) but seems 
to be new for the Vlasov-Darwin system: we consider classical solutions of the Cauchy 
problem and show that these exist for all times if the initial data are chosen sufficiently 
small. The precise statement of our result is contained in the next section, where we also 
formulate three propositions which are used to prove the theorem. Sections [3H6] are devoted 
to proofs. In the final Section[7]we take a look at spherically symmetric solutions. First it is 
shown that any symmetry of the initial datum with respect to an orthogonal transformation 
is preserved for all times. This allows the conclusion that in case of spherical symmetry 
the VD system reduces to the well known relativistic Vlasov-Poisson system. So in this 
case the solutions are global-in-time as well [5J. 




(1.3) 



E = E L + E T , VxE L = 0, V-E T = 0. 



(1.4) 



d t E L -VxB = -j, V-E L = p 
d t B + V xE T = 0, V-B = 0. 



(1.5) 
(1.6) 
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2 Results 



Before presenting the main result of the present paper (which is formulated as the following 
theorem) we fix some notation. Let Rq,Pq>0 be fixed throughout the paper. For r > let 
B r := -B r (0) = {x £ M 3 : \ x\ < r}. Furthermore we specify the set where the initial data are 
taken from: Let C 2 (W 1 ) denote the space of twice continously differentiable functions on 
M. n with compact support and 

V := {f G C 2 (M 6 ) : /> 0, ||/][oo < 1, \\VfU < l,supp / C B Ro (0) x B Po (0)}. 

The Legesgue space of square integrable functions is denoted by L 2 (M 3 ) and P: L 2 (]R 3 ) — ► 
L 2 (K 3 ) is the projection on the divergence free part, which is discussed in the Appendix. 

If I is an intervall and g: IxW 1 — >W n , we denote the quantity sup x£K n\g(t,x)\ by 
||<?(*)||ooj and if K cR n , the expression ||g(t)||oo,_ft: means sup^g^ \g(t,x)\. 

Theorem 2.1 There exists 5>0 such that the classical solution of the VD system with 
initial datum fo in T> satsifying ||/o||oo<^ exists globally in time. 

For the proof of this result the reformulation of the field equations of the VD system 
in terms of potentials $ and A given in [8] is used. Let 

A$ = p, lim $(i,x) = 0, (2.7) 

|x|^oo 

AA = -F(j), lim A(t,x) = 0. (2.8) 

\x\^oo 

Then the components of the electromagnetic field are 

E L = V®, B = VxA, E T = -d t A, 

cf. Lemma 2.3 in [8]. 

The proof of the theorem is given in Section [6l Sections O [U [5] contain the proofs of 
preliminary results, which are formulated in Propositions 12.31 12.41 and 12.51 A prominent 
role in the following is played by a certain decay condition. 

Definition 2.2 A classical solution f : [0,T[xR 6 ^lR of the VD system is said to satisfy 
a free streaming condition with parameter a on an interval [0,a] if 

\\E T (t)\\ 00 + ||£ L (t)||oo + ||5(t)|| 00 <a(l + t)- 3 / 2 , 
\\VE T (t)\\ 00 + IIV^^IIoo + UVB^IIoo^aa + t)- 5 / 2 

for all t G [0,a]. 
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As for the RVM there is a continuation criterion for solutions of the VD system, 
which says that solutions may be continued as long as the momentum support {pSl 3 : 
3x,t: f(t,x,p) 7^0} remains bounded [8j. Using this criterion one can show easily that a 
solution which satisfies a condition as the one above on its maximal interval of existence 
is indeed a global one. 

We now discuss the main idea of the proof. The task will be to show that the free 
streaming condition implies decay of the source terms p and j as well as decay of the fields 
El,Et,B. This will be done in the following two propositions. 

Proposition 2.3 There exist a, Ci(Rq,Po) > such that for every solution of the VD 
system with /(0)gP which satisfies a free streaming condition on an interval [0,a] with 
parameter a it holds 

iipWHoo + ii j moo < Ciitr 3 , iim*)Hoo + \\9 x m\\oo < c x . 

Proposition 2.4 There exist a, C2(Ro,Po) > such that for every solution of the VD 
system with /(0)eP which satisfies a free streaming condition on an interval [0,a] with 
parameter a it holds for t£ [l,a] 

||^T(i)lloo + ||^(i)lloo + ||^(t)||oo<C 2 t- 9 / 5 , 

||V£ T (i)||oo + \\VE L (t)\\ oa + 1| VtfWHoo < C 2 t- 8/3 

These estimates provide the main ingredient needed for the bootstrap argument in the 
proof of the theorem. It follows that a solution satisfying a free streaming condition decays 
asymptotically even faster. This is an important point for the global existence argument. 
To start this bootstrapping we need a further tool which is given in the next proposition: 
if the initial datum is chosen sufficiently small, then the fields remain small for some time. 
This may be interpreted as a statement about continous dependence on initial data. 

Proposition 2.5 Let e,T>0 be given. Then there exists 5>0 such that any classical 
solution f of the VD systems with /(0)e£> and \\f(0)\\oo < 5 exists at least up to time T 
and is such that 

\\EL(t)\\oo + \\E T (t)\\oo + \\B(t)\\ 00 + \\VE L (t)\\ 00 + ||V£ r (t)||oo + ||V5(t)||oo <e. 



4 



3 Decay of the source terms 



For proving Proposition 12.31 we need the following 

Lemma 3.1 Lett>0, £eC 2 ([0,t]), = £(£) = () and let 

|^)l<ci( S ) + c 2 ( S )|C( S )| + C3 ( S )|e(s)| 
for s£ [0,t], where ci,C2,C3>0 are continous and C3 is monotonically decreasing. Then 



Proof. Define z(s) := J*|C(r)|dr, so that |f (s)| < z(s) = -\£(s)\,z(t) = z(t) = 0. 



Obviously 



z{s) 



'({a)da 



dr 



<J s J \i(<r)\dadr, 



so that by our assumptions 

ft ft ft ft ft ft 

z(s) < ci(a)dadr+ / / C2(cr)z(a)dadr — / / c^(a)z(a)dadr 

Js Jt Js Jt Js Jt 

= : h + I 2 +h. 
Changing the order of integration it follows that 



ft rcr ft 

I I ci(a)drda< / ac\(a)da 

J s J s J s 



h 



as well as 



t per 




C2{o-)z(a)drda < / ac2(c)z(a)da. 



The integral ^3 can be estimated in the following way: 

L 3-J s c s( r ) (- J z(o-)da^j dr = j c 3 (T)z(T)dr, 

where the monotonicity of C3 was used in the first step and the relation z(t) = in the 
second. Hence the function z satisfies the integral inequality 

z(s)< / aci(a)da+ / [ac 2 (a) + 03(a)] z(a)da, 

J s J s 



so that by Gronwall's Lemma 

<*(s) < (J aci(a)daj e X>W CT )+ c 3(*. 

□ 

Proof of Proposition \2.3l 

The proof presented here is an adaptation of the corresponding argument for the Vlasov- 
Poisson system (cf . [9] ) . To get decay of p a change of variables is performed in the integral 
defining it. The transformation determinant appearing can be shown to decay fast enough. 

Let / be a classical solution of the VD system with f(0) = f°£V and denote by 
(X(s,t,x,p), P(s,t,x,p)) the corresponding solution of the characteristic system 

X(s,t,x,p) = v(P(s,t,x,p)), 

P(s,t,x,p) = E(X(s,t,x,p),s) + v(P(s,t,x,p)) x B(X(s,t,x,p),s) 

with initial condition X(t,t,x,p) = x,P(t,t,x,p) =p. 
Then we have 

f(t,x,p) = f°(X(0,t,x,p),P(0,t,x,p)). 
If \p\ <Po, then by the free streaming condition 

\P(t,Q,x,p)\<P + [ (||^(s)||oo + ||B(s)|| oo )ds<P + a f (l + s)- 3/2 ds<P + 2a. 
io Jo 

Hence for small enough a we may conclude that \p\ >-Po + l implies f(t,x,p) = 0. 
Define 

£0) :=d p X(s,t,x,p)-(s-t)Dv(p). 
We have = 0, and using the characteristic system one obtains 

i(s) = Dv{P{s))d p P{s) - Dv(p). 

So £(t) = 0. Differentiating once more we get 

Z(s) = D 2 v(P(s))P(s)d p P(s) + Dv(P(s))d p P(s). 

As is easily checked, Dv(p) and D 2 v{p) are bounded independently of p, so that 

m\<c(\P(s)\\d p p(s)\ + \d p p(s) 

and therefore by the characteristic system 

ie*(^)l < 0(1^^(^)1(11^(^)1100 -h 11^(^)1100)+ (11^^(^)1100 + 11 v^C^)!!^)!^^^)!) . 
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Resubstituting we have 
and 



\d p X(s)\<\(i(s) + (s-t)Dv(p)\ 



d p P(s) = (DviPis)))- 1 (i(s) + D p v). 
Assuming \p\ <-Po + l> we can estimate 

\d p P(s)\<C(P )(\i(s)\ + l 



Using the free streaming condition, we finally obtain the following second order differential 
inequality for £: 

||( S )|<C(Po)a{(l + sr 3/2 + (t^^ 
By Lemma 13.11 

\as)\<C(P )a(t-s)e c ^ a , 

where we possibly have to adjust the constant C{Pq). In terms the characteristic variables 
this means for a chosen sufficiently small that 

\d p X(0,t,x,p) + tDv(p)\ <et, 

where e > is prescribed such that 

\p®p 



, .„ +e v / l + |p| 2 </3<l for \p\ <P + 1. 
l + \p\ z 

Here p0p denotes the matrix whose (i,j')-entry is PiPj. This implies 



d p X(0,t,x,p) + 



By direct computation 



1 

Dv{p) 



<et + 



p®p 



I-tDv(p) 



hence 



d p x(o,t,x, P ) + 



V 1 +\p\ : 



<et + 



t \p®p\ t 

< , J - 



So the linear map d p X(0,t,x,p) is invertible and in conclusion the transformation 
\E': Bp 0+ i -^M 3 ,wi-^ X(0,t,x,p) is a local diffeomorphism. It is even a diffeomorphism onto 
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its image, since it is ono-to-one as well: Let p,p£ -Bp +i be given and p T :=rp+(l — r)p. 
Then 



|*(p)-*(p)| 



<9 p X(0, t,x,p T )(p-p)dr 



d p X(0,t,x,p T ) + 



V 1 + \pt 



(p-p)dr- 



f 1 1 
> t\p-p\ / -7= 

Jo + 



>C(l-(3)t\p-Pl 



-dr-(3t\p-p\ 



Vi + brP 



(p — p)dT 



where C depends only on Po- 

Denote the open range of ^ by f/ and let <3?: ?7^-Bp 0+ i be its inverse. Calculation 
gives 



p(t,x) = / f(t,x,p)dp 

•> b p +i 



-J 
= / 



/°(tf(p),P((M,z,p))dp 



<£>(*(Bp ()+ i)) 

f°(w,P(0,t,x,$(w)))\detD<5>(w)\dw. 



'*(Bp 0+ i) 

For the functional determinant showing up here we have by our previous calculations 

1 



| detD$H| = | det[.Dtf 1 



|detD*($H)| 



(3.9) 



Note that 



det£>*0) = det d p X(0,i,x,p) + 



(l + |p|2)3/2 



det \A + |p| 2 - 



d p x(o,t,x, P ) + 



But for the first matrix in the argument of the determinant we have 



_d p x(o,t,x,p) + 



sfi+\p\ 



-I 



</?<!, 



8 



so the absolute value of the determinant is bounded from below by a positive constant Cp. 
Returning to (|3.9|) . it is seen that 

|detD*(«;)|<^^, 

resulting in 



t 3 

In addition 

C* 

\j(t,x)\<p(t,x)<—. 
For the bounds to be obtained for d x p and d x j note that 
\d xP (t,x)\ <C(P + l) 3 ||^./(t)||oo, 

\d x f(t,x, P )\ < c(\d x x(o,t,x,p)\ + \d x p(o,t,x, P )\). 

Next let £(s) := d x X(s,t,x,v) - id such that £(s) = Dv(P(s))d x P(s). Obviously f (i) = £(t) ■■ 
0. Differentiating further one has 

|£(s)| = \D 2 v(P(s))P(s,t,x,p)d x P(s,t,x,p)+Dv(P(s))d x P(s,t,x,p)\ 
<aC((l + s)- 3 / 2 \d x P(s,t,x,p)\ + (l + s)- 5 / 2 \d x X(s,t,x,p){ 



where again the decay of the fields due to (almost) free streaming was employed. By 
definition 

\d x X( 8 )<\S( 8 )\ + l and \d x P(s)\<C\£(s)\ 
and we may assume \p\ <Po + l to discover the relation 

||( S )|<Ca{(l + S )- 5 / 2 |^)| + (l + S )- 3/2 |^( S )| + a(l + S )- 5 / 2 }, 

which by Lemma 13. II implies 

\as)\<Caj\l + a)-^ 2 dae Ca f^ 1+ ^- 3/2 ^<2Cae 2Ca ^ 
An easy application of GronwalFs Lemma shows that |£(s)| is bounded too, which means 

(\d x x(o,t,x,p)\ + \d x p(o,t,x, P )\) < c 

and all claims are proved. □ 
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4 Decay of the fields 



Proof of Proposition \2.4\ First let the constants a and C\ be as given by Proposition | 
We want to get sufficiently good decay rate estimates for the fields from the field equations, 
the decay of the source terms, and the free streaming condition. The field Et is treated 
first. We have 

AE T = -d t (AA) = d t (¥j) =¥(d t j), 

where the last equation, i.e. the commutativity of P and dt, is read off directly from the 
Fourier representation of the projection operator, compare (|4.10p or the Appendix. 
The Vlasov equation, Eq. (jl.ip . then implies 

d t j(t,x) = J d t f(t,x,p)v(p)dp = - J (v(p),V x f)v(p)dp- J v(p)®K(t,x,p)V p fdp. 

Here K(t,x,p) = E(t,x) + v{p) x B(t,x). Integration by parts in the last term finally leads 
to 

div x [f(t,x,p)v(p)®v(p)]dp+ / - f(t,x,p)K(t,x,p)dp 

J vi+P 

=: G 1 (t,x) + G 2 (t,x), 

where the divergence appearing is to be understood row-wise. 

Writing Et = E T + E T , where the components of the r.h.s. are solutions of AE T = 
F(Gi) and AE T = W(G2) respectively, we treat each of them seperately. Recall the Fourier 
representation of the projection operator P: 

¥F(x) = J e ^ k?I ~^ k F{k)dk, (4.10) 



F(k) = (2tt)- 3 J e-' lkx F(x)dx (4.11) 
(compare the Appendix) . The solution of the Poisson equation may then be expressed as 

EUt,x) = -je^ lkl2l -^ k G l (k)dk, 1 = 1,2. 
Introducing M = (Mi,M2,M^) = J f(t,x,p)v{p)®v{p)dp, it comes 
E T (t,x) = Je^ lk2lI l - k ® k G 1 (t,k)dk 

= Ej^ xlJ %^^M j{ t,k)dk 
= iJ2 I d kx mj(k)Mj{t,k)dk, 
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where nij{k) = ^ ^ffcp^ Tfci * s a f unc tion homogeneous of degree -1. The theory of pseudo 
differential operators (compare (45) and the proof of (47) in [7]) permits us to estimate as 
follows: 

WE^W^KCiWMW^ + WMWp), l<p<3, (4.12) 
||V x ^(t)||oo < C(7 3/p '- 2 ||VM|| p + log( 7 - 1 )||M|| 00 + ||M|| g ), (4.13) 

3 <p < oo, 1 < q< oo. Here the parameter 7 is restricted to the interval ]0,1]. Since 
(t) Hoc < ||p(t)||oo < Ct~ 3 , the first equation implies setting p = 5/2: 

ll^(i)||oo<Ct- 9 / 5 



Using estimates from the proof of Proposition I2T31 and HVM^Hqq < ||Vp(i)[|oo <C it follows 
that i 

||VM(t)|| p < ||VM||^r||VM||} /p <C||VM||} /p <Ct 3 / p , 
so that setting 7 = t~ 3 we get for t > 1: 



(here the choice p = 36,g = 17/9 was made). 
Now consider AE^ = F(G 2 ). We have 

E 2 T {t,x) = j e ikx ^ 2l ^® k G 2 (t,k)dk = J e lkx m. 2 G 2 (t,k)dk, 

d 3 E 2 T (t,x) = i f e lkx lkl2l ~^ k ^G 2 (t,k)dk = J e tkx m^G 2 (t,k)dk. 

The symbols m a showing up here are homogeneous of degree a. A simple adaptation of 
the proof of (45) in [7] shows that 

||^(t)||oo<C(||G 2 (t)|| 0O + ||G 2 (t)|| p ), l<p<3/2, (4.14) 

and as before 

IIS^-E^Wlloo < C(l|G r 2(*)Hoo + II^WHp), 1<P<3. (4.15) 
Since \G 2 (t,x)\ < WKtyWoopfyjx), we have 

||G 2 (i)|| 5/4 <C7(l+i)- 3 / 2 ||p(t)||V 5 ||p(t)||J /5 <r 21 A0 

and 

||G 2 (t)|| 5/2 <C(l + t)- 3 / 2 ||p(t)||^ 5 ||p(t)|| 2 / 5 <t- 33 / 10 . 



11 



Altogether this implies 

\\E T {t)\\oo<Ct-V\ 
||V^ T (t)||oo <Ct~ & l\ t>l. 

Now we come to the other fields. The longitudinal part El of the electric field is treated 
exactly as in the case of the Vlasov-Poisson system (cf. [S]): 

\\E L (t)\\oo <Ct~\ 

\\d x E L (t)\\oc <Ct- 3 logt. 

The bounds for the magnetic field B = V x A field are obtained in a way analogous to 
the procedure used so far. First we have a representation 

VA(t,x) = J e^m-iWj&fydk 

with m_i homogeneous of degree —1. Therefore 

\\B(t)\\oo<\\m\\oo+\\mh/2<ct- 9 / 5 . 

In analogy to our treatment of Ej- we find 

||V a; S(t)|| 00 <c(7 3/p '- 2 ||Vi|| ? ,+log(7- 1 )||i||oo + ||i|| ? ) (4.16) 
and the proof may be completed as shown before. □ 

5 Continous dependence 

In this section we denote by C a constant depending only on Rq,Pq, which may change 
from line to line. For the proof we collect some facts first. 

Let (f,EL,ET,B) be a solution of the VD system on some time interval [0,T[ with 
/o = /(0)€ V. Define 

Q{t):=sup{\p\:3x,0<s<t: f(t,x,p)^0}. 
Then we have the following 

Lemma 5.1 Let f be a solution with f(0) £T>. Then there holds 

\\p(t)\U/3 + \\j(t)\U/3<C\\fo\U 

PWIIoo^CH/ollQW 1 / 3 , 
\\VA(jb)\\oo + \\V*V)\\oo < C||/ ||ooQ(t) 5/3 . 
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To estimate the field Et only a local result is available. 

Lemma 5.2 Let f : [0,T] xl 6 ^l be a solution with f(0)eT>. Then 

||E T (t)|| 00iBflo+T <C Bo+T (l + ||p(i)|| 3 )(||F(t)|| 6/5 + || J P(i)|| 2 ) 
where F = F\ + F2 with 



F x {t,x) = J div( x )(f(t,x,p)v(p)®v(p))dp, 

F 2 (t,x) = [ I ~ V ^l^ P) f(t,x,p)(E L (t,x)+v(p)xB(t,x))dp. 

J V 1 + \p\ 

Detailed proofs of these Lemmas are given in Pallard's paper [8], where also the following 
theorem is proved. 

Theorem 5.3 Lei f £C 2 (R 6 ). Then there exists T* > and a unique solution 
(f,EL,ET,B) to the Vlasow-Darwin system with /(0) = /o satisfying 

f G CHlO^fxIR 3 xK 3 ), 
E L ,B G (^([Cr-fxM 3 ), 
E T ,V x E T eC([0,T*[xR 3 ), 

and such that for any t£ [0,T*[ the distribution function f(t,.) is compactly supported. 

By inspection of the constants in the proof one finds that a strict lower bound for T* 
is given by T 1 := (C||/o||^o) _1 with a constant C independent of fo, i.e T* > T' . In addition 
one has 

Q(t)<C(R ,P ) (0<t<T'). (5.17) 
One last ingredient for the proof of Proposition 12.51 is contained in the following 

Lemma 5.4 Let f be a solution and T' be defined as above. Then 

\\V x , P f(t)\\oo<C(Ro,P ) (0<t<T'). 

The proof can again be found in [8J. From the Lemma we immediately deduce the 
bounds 

||V^(t)||oo + 1| Vj(t)||oo < C7, 0<t<T'. 
With these tools at hand we are now ready for the 
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Proof of Proposition \2.5[ 

Let e,T>0 be given. From the above facts one finds immediately that the solution 
interval can be made as long as we wish and that ||-Ex(t)||oo and ||-B(i)||oo can be made as 
small as necessary by choosing 5 sufficiently small. 

It is standard to obtain a bound for llV^-K^^Hoo (see [9J), and for ||V-B(i)||oo we can 
use Eq. (|4.16p from Section [U By finite propagation speed and since ||Vj(t)||oo remains 
bounded for ie[0,T"], we can choose the parameter 7 on the right-hand side of the in- 
equality properly to get the result. 

We still have to get control over ||£r||oo and ||VEt||oo- From (|4.12j) and (|4.14p we have 
the estimate 

pT(t)|| o<C(||p(t)|U + ||p(t)|| 2 + ||G 2 (t)|| QO + ||G 2 (t)|| 5/4 ) J 
where the notation introduced in Section 2] is used again. Now 

I — v(p)®v{p) 



\G 2 (t,x)\ < 



-f(t,x,p)K(t,x,p)dp 



<C f f(t,x,p)\K(t,x,p)\dp 
J\ P \<Q(t) 

<CQ(tf\\f \U\K(t)\\oo,B RQ+T XB Ro+T ^)- 



and since we have bounds for ||^(i)||oo,_B flo+T by Lemma [5T2l we get 1 1 -Z?^- 1 (t ) 1 1 00 < e for ||/( 



0||oo 



chosen sufficiently small in view of (|5.17p . 

To estimate ||V£Jt(*)|| we note as in Section HJ 

Hv^^c*)!!^ < c r-r 3 /^'~ 2 ||VAi r || i:> +iog(^- :L )||iv^|| 00 + ||a^|| 2 + 11^2(^)110.3 + ||C 



with 0<7<1,3<j>< oo. So again each term can be made as small as wished and the proof 
is complete. □ 



6 Proof of the theorem 



We start by choosing a constant To > such that for t > To it holds 

C 2 t- g/5 <a(l + t)- 3 / 2 and C 2 t- 8/3 < a(l + t)~ 5 / 2 , (6.18) 

where a and C 2 are the constants given by Proposition 12.41 Proposition 12.51 says that 
there exists 5 > such that a solution of the Vlasow-Darwin system with initial fo€T> and 
1 1 /o I |oo < S satisfies 

\\E L (t)\\ 00 + \\E T (t)\\ 00 + \\B(t)\\ 00 + ||V£k(t)|| 00 + ||V#r(*)||oQ + ||Vfl(t)||oo 

<(1 + T )- 5 / 2 , 
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for t belonging to [0,To + l]. Moreover, it may be assumed that the maximal interval of 
existence 1= [0,T max [ is strictly larger than [0,To + l], i.e. To + l<T max . 

If / is a solution as above, then by continuity / satisfies a free streaming condition with 
parameter a on an interval [0,T*] with Tq<T* < T max and T* may be chosen maximal with 
these properties. Because of equation (|6.18p we may now conclude with Proposition 12.41 
that 

||^L(t)||oo + ||^T(t)||oo + ||5(t)||oo<a(l + t)- 3/2 , 

HV^lWHoo + l|V£ T (t)||oc + ||V£(£)||oo < a(l + t)- 5/2 
for all t£l. But this implies T max = oo and the solution is global. □ 

7 Spherically symmetric initial data 

In case the initial datum f° is spherically symmetric, which in the present situation by 
definition means 

f°(Qx,Q P ) = f°(x, P ) \/x, P eR 3 ,Qe 0(3), 

the Vlasov-Darwin system reduces to the relativistic Vlasov-Poisson system, as is shown 
in the following. First we show that spherical symmetry is preserved. 

Lemma 7.1 Let f : [0,T[xM 6 ^]R be a classical solution of the Vlasov-Darwin system and 
let /(0) be spherically symmetric. Then f(t) is spherically symmetric for all 0<t<T. 

Actually we will see that any invariance of the initial datum with respect to an orthog- 
onal transformation is preserved for all times, which implies at the same time that e.g. 
cylindrical symmetry or reflectional symmetries are preserved as well. 

Proof of the Lemma. 

Let QgO(3) be given and set f(t,x,p):=f(t,Qx,Qp). It suffices to show that / solves 
the Vlasov-Darwin system. One finds 



p(t,x) := J f(t,x,p)dp = p(t,Qx), 

j(t,x) := / f(t,x,p)v(p)dp = Q~ 1 j(t,Qx). 



For the potentials <j> and A one therefore has 



$(t,x) = &(t,Qx), 
A(t,x) = Q^A(t,Qx), 
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as can be seen easily from the Fourier representation of the projection operator P, compare 
P~TU|h This implies 

E L (t,x) :=V$(t,x) = Q- 1 E L (t,Qx), 
E T (t,x) := -d t A(t,x) = Q- 1 E T ( y t,Qx). 

We set B := V x A. Then by Lemma 2.3 in [Hj and since 

d t p(t,x)+V-j(t,x) = d t p( y t,Qx)+V-j(t,Qx) = 

the quantities (El,Et,B) solve the field equation part of the Vlasov-Darwin system. 

We have to show that the transport equation (jl.ip holds. Consider the term 
p x B(t,x) = px (V x A(t,x)). By well known vector identities we can write 

pxB(t,x) = S7{A(t,x)p)- (p-V)A(t,x)) 

p. 



DA(t,x)j -DA(t,x) 

Here D denotes the total derivative w.r.t. x. Now DA(t,x) = Q- 1 DA(t,Qx)Q and therefore 

Q(pxB(t,x)) = [(DAfaQx))* - DA(t,Qx)]Qp 
= QpxB(t,Qx). 

The last equality holds because the forgoing applies equally well to A as to A. This finally 
leads to t£TQ). □ 

So we have seen that spherical symmetry is preserved for all times. This includes that 
for te [0,T[,Qe 0(3) the following identities hold 

p{t,Qx) = p{t,x), 
$(t,Qx) = $(t,x), 
j(t,Qx) = Qj(t,x), 
A(t,Qx) = QA(t,x). 

Lemma 7.2 The vector field j is radial. 

Proof. In the following the dependence of j on t is suppressed. Let x E M 3 \ {0} be given 
and choose a positive orthonormal basis (61,^2,^3) with b± = Let Q\,Q2 be orthogonal 
transformations with matrices 

Mi 





0) 




(-10 




• 


-1 


, M 2 =\ 


1 


• 


\ 


U 




V 
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w.r.t. to the basis chosen. Let j(x) = Y^,j a j^j- Then 

Qdi x ) = -a.\b\ - a 2 b 2 + 0363, Q2j{x) = -0161 + 0262 - 0363. 
But since Qij(x) = j(Q\x) = j{— x) =j{Q2x) = Q2j{x) it follows that 02 = 03 = 0. □ 

Lemma 7.3 It holds Pj = 0. 

Proof. First note, that V ■ j(Qx) =V ■ j(x) for all Q G 0(3). Recall the definition of P: 

Fj(x)=j(x)+V^(x) 

where 



47T 7 p — J/ 1 

But since the source term V • j has rotational symmetry, the forgoing simplifies to 

f r x 
W(x) = - / s 2 (V-j)(s)ds-3, r = |x|. 

Jo r 
The integral in the last expression can be transformed to 

r s 2 (v.j)( S )d S =^ j v-jdv 

4tt J dBr 
= j(x)nr 2 . 

So 

x 

aj) = + V^fx) =j(x) — j(x)n— = 0. 

r 



□ 



Lemma l7.3l implies that also A(t) = for t 6 [0,T[. This immediately leads to Et = B = 0, 
so that the following proposition is proved. 

Proposition 7.4 In case of a spherically symmetric initial datum f° the Vlasov- Darwin 
system reduces to the (spherically symmetric) relativistic Vlasov-Poisson system (with re- 
pelling forces). Hence in this case the solution is global. 

The proof of the last statement is given in [5] . 
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Appendix 



We start with some remarks about the projection operator P: L (Mr) L (Mr) , which is 
defined as follows: For F£C^(R 3 ;R 3 ) one sets 

(FF)(x) = F(x) + V$(x), where 



4-7T J \x — y\ 

Since — = X7 ■ F we cleary obtain V-P-F = 0. Applying the Fourier transform to these 
relations it follows 

VF(0 = F(0 + iHm, 

hence 

MO={l~^f)HO- (7-19) 

Therefore |P.F(£)| <C\F(£)\, so that by the Plancherel-Theorem P extends to a continous 
operator on L 2 (1R 3 ) characterized by (17.190 . 

We conclude with some remarks about the pseudo differential operators used in Section 
[H Such an operator is of the form 

Au(x) = j^J MOHO^M, (7-20) 

where n is the dimension of the underlying space 1", the function Aq is called the symbol 
of the operator and is chosen from a suitable set of functions, and u is the Fourier transform 
of u. We can restrict ourselves to the case that u belongs to the Schwarz space 5(M n ) of 
rapidly decreasing functions. It is shown e.g. in [4] that an operator of the form (|7.20p 
with a symbol Aq homogeneous of degree a>— n, i.e. Ao(tt;) = t a Ao(£) for i>0,(Gl", 
has a representation as an integral operator of the form 

Au(x)= / a Q (x-y)u(y)dy, 



where ao is a function homogeneous of degree —a — n. For such a (smooth) function one 
clearly has 

\a (y)\<C\y\- a - n 

and this is all one needs to know for the estimates presented here and in [7j. 

Acknowledgement. I would like to thank Prof. G. Rein for bringing me in contact 
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